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The discontinuous nature of chromospheric activity 
evolution 
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Abstract Chromospheric activity has been thought to 
decay smoothly with time and, hence, to be a viable 
age indicator. Measurements in solar type stars in open 
clusters seem to point to a different conclusion: chromo- 
spheric activity undergoes a fast transition from Hyades 
level to that of the Sun after about 1 Gyr of main- 
sequence lifetime and any decaying trend before or af- 
ter this transition must be much less significant than 
the short term variations. 
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Introduction 

In the atmosphere of the Sun and of the stars similar 
to it, several processes take place, which make the tem- 
perature of the atmosphere higher than it would be if 
radiative equilibrium held. Such non-radiative heating 
mechanisms are powered by convection and magnetic 
field. In the outer layers of these stars, the temper- 
ature is increasing towards the surface, and the main 
cooling mechanism is radiative loss through strong res- 
onance lines, such as C a II H and K, Mg II h and k, 
Hq,. (see e.g. Hall 20081 for a review). This region is 
referred to as chromosphere, the radiation emitted in 
it is called chromospheric activity (hereafter CA), and 
the strength of the emission cores in some of the afore- 
mentioned resonance lines are a good CA indicator. 

Solar-type stars loose angular momentum due to 
magnetised stellar winds, and CA is believed to decay 
with time as a consequence. This trend is superim- 
posed to short term CA variations, such as those that 
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caused the Maundcr-minimum in the Sun, activity cy- 
cles equivalent to the 11-year long solar o ne, and, in the 
most active stars , rotational modulation rtWilsonlll963 ; 



Skumanichlll972l : iNoves et alJll984t ISimon et al 



19851) 



Due to its decay, CA is a potential age indica- 
tor and several efforts have been undertaken in or- 
der to calibrate it (see 





Lachaume et al.l Il999t iMamaiek fc Hillenbrand 



e^ 



Soderblom et al 





1991 




2008 



and references therein). But results coming from the 
analysis of solar type stars in open clusters are chal- 
lenging the view of a smooth, decreasing trend modu- 
lating CA cycles. We do not question the main picture 
described above on CA decay, but it appears that this 
occurs mainly in one major event, in which it spans the 
Vaughan-Prcston gap, and any evolution before and af- 
ter such event is less important than the variation due 
to activity cycles. 

Rotation i s probably better correlated with stella r 
ages (see e.g. lBarnesll2009l : ICardini fc Cassatellall2007l ). 
since it is not affected by any cyclic variation. How- 
ever, rotation periods measurements are not possible 
for inactive stars, for which we can only have an esti- 
mation of the projected rotation velocity, which is not 
sufficient to to draw any firm conclusion on the nature 
of the angular momentum evolution. 

The matter deserves close attention, because chro- 
mospheric ages have been used in important studies on 
chemical enrichment and star form ation in the Galac- 
tic disk jRocha-Pinto et ah 2000alfbh, and the a ge of 
stars hosting planetary systems (jSaffe et al.ll2005l) . We 
should reach a definitive conclusion on whether CA is 
or not a good age indicator, and this can be better 
achieved by observing more solar-type stars in open 
clusters and by determining more precisely their ages. 

Intense investigation has long been conducted on the 
topic, especially by using the large a mount of dat a col- 
lected in Mount Wilson ca mpaign rtWilsoiJll96"8h and 
in the planet-search surveys (Isaacson 20091) . However, 
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only the possibility of observing at high signal to noise 
and high resolution solar type stars in old and distant 
open clusters, provided by the 8- and 10-m class tele- 
scopes, could make us progress in this field, and it is 
the key to address still unanswered questions. 

In the present work, I will report on the results that 
led us to reconsider the viability of CA as age indica- 
tor. It is also shown that data used to calibrate chro- 
mosphcric ages do not contradict our conclusions. 



1 Data sample 

The complete sample on which our conclusions are 
based consists of 40 main-sequence stars in 7 open 
clusters and the Sun. The nature of single, dwarf 
members of the parent clusters were established, for 
the 40 target stars, in published p hotometric, proper- 
motion and radial-velocity studies (IHerzog et al.lll975 : 



Latham et al. 1992; Nordstroem et al. 1996; Mcibom et al 



20021 iMermilliod k Mavorl [l99a iTwarog et ailll993l) . 
and, only for IC 4756 and N GC 5822, from ou r own 



radial- velocity measurements ( Pace et al. 2010l) . The 



data about 15 stars in the Hyades cluster consist of 
spectra taken with HIRES at Keck, kindly made avail- 
able to us by D. Paulson, A. Hatzes and P. Cochran. 
The rest of the data consists of spectra taken with 
UVES at VLT in two different runs. In the first (ESO 
run 66D-0457, P.I. L. Pasquini) the targets were: 7 in 
Praesepe, 2 in NGC 3680, 5 in IC 4651, and 6 in M 67. 
In the second run (ESO run 73D-0655, P.I. G. Pace) 2 
stars were observed in NGC 5822 and 3 in IC 4756. 

Our UVES spectra have a resolution of R«100 000 in 
the spectral range from 4 800 to 6 800 A, and R«60 000 
in the spectral range 3 200 to 4 600 A. After summing 
the spectra of the same star, we achieved S/N ratios 
per pixel ranging from about 50 to about 150. Due to 
the higher apparent brightness of its stars, the quality 
of Praesepe spectra is remarkably high. The Hyades 
spectra we used have a resolution of R~60 000 and a 
signal-to-noise ratio per pixel ranging from 100 to 200, 
and from 20 to 30 at the core of the Ca II K line. 



2 Data analysis 

The procedure we adopted to measur e chromospheric 
activ i ty is described in detail elsewhere ( Pace k Pasquini^ 
20041 Pace et al]l2009t ) . Here I present a complete sum- 
mary. 

As an indicator of chromospheric activity we used 
F' K , i.e. the energy flux of the Ca II K line emitted per 
unit surface in the chromosphere. We also computed 
the value normalised to bolometric emission: 



\ogR' K = \og[F' K /aT^ }j 

The spectra were normalised by dividing their in- 
tensity by the counts at the pseudo continuum point 
3950.5-A. 

The normalised flux was integrated over a 1-A wide 
region centred on the Ca II K line. This region co- 
incides with the chromospheric emission peak. The 
result of this integration is the 1-A K index and in- 
cludes a contribution from the photosphere. The stars 
in NGC 3680, IC 4651, M 67, IC 4756 and NGC 5822 
showed the interstellar K absorption line which affected 
the chromospheric K-line feature. As for NGC 3680, 
IC 4651, and M 67, we could evaluate the contribution 
of the interstellar absorption, by observing hot stars in 
each of these clusters, and correct for it. This was not 
possible, instead, for stars in IC 4756 and NGC 5822. 
To measure the 1-A K index in these stars, we inte- 
grated the normalised flux only in the uncontaminated 
-part of their profile. The Praesepe stars are unaffected 
May interstellar absorption thus allowing us to calculate 
a ratio between the flux measured over the 1-A and 
that measured over the portion of the feature which is 
unaffected by interstellar line in all stars. The initial 
measures for the affected stars were then multiplied by 
this factor to give a final corrected 1-A K index. The 
errors involved in the measurement of the 1-A K index 
were evaluated to be within 6%. 

For the Sun, we used the 1-A K index measurements 



made bv lWhite k Livingston! (|1981I ). They monitored 
solar chromospheric activity from the first minimum to 
the maximum of the 21 st solar-activity cycle. 

Subtraction of the photospheric contribution to the 
1-A K index was made as follows. For the Sun, we 
computed the photospheric correction using a solar- 
photosphere synthetic spectrum (courtesy of P. Bonifa- 
cio). For the other stars, the photospheric contribution 
was computed by scaling the solar photospheric con- 
tribution by a factor that depends on the stellar pa- 
rameters. These scaling factors were co mputed using 
the spectral synthesis code of MOOG dSnedenl Il973l 
versio n 2002), and Kurucz's grid of models ( Kurucj 



1993). Stellar parameters, namely temperature, grav- 



ity, microturbulence and metallicity. were known from 
our chemical analyses in iPace et al. 1 d2008l l2010h and, 
for Hyades stars, from lPaulson et al] (|2003f ). 

In order to transform the 1-A K index, which is an 
equivalent width, into an intrinsic flux, namely F' K , we 
multiplied the former by the flux at the stellar surface 
of the pseudo-continuum radiation at 3950.5-A, which 
was obtained from the stellar tempe r ature s using re- 
lation ships published in ISousa et al.l (|2008l ); iPasquini 
( 19851) ; and lPasauini et al.l (|l988l) . 

The employment of stellar temperatures from spec- 
troscopic analysis of iron lines instead of published 
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colours avoids the use of uncertain reddening estima- 
tions. 

The results of the reviewed analysis of the old sample 
and of the analy sis of the new sample, published in 



the observed flux of the inactive stars. In order to pro- 
vide a physical exp l anatio n for the occurrence of the VP 



Pace et al.l (|2009[ ). were used to produce the diagram 



of CA as a function of temperature shown in Figure [TJ 



3 Discussion. 

It was already pointed out back in the eighties that the 
distri bution of CA is markedly bimodal (jVaughan fe Presto 
1980). Namely very few stars are less active than the 



gap, iDurnev et al.l (|198ll) proposed a transition from a 
complex to a simpler magnetic-field morphology which 
occurs at the time when the rotation decreases enough 
to rea ch a threshold value. More recently, iBarnes 
( 20031) detected two sequences in the pcriod-versus- 
colour diagram of open clusters, and he associated them 
with two different rotation morphologies, i ntertw ined 
with stellar magnetic fields. iBohm-Vitensd (|2007t ) sug- 
ested a change of dynamo mechanism to explain the 



ges 



Hyades and significantly more active than the Sun. 
This underpopulated range of values i s usually referred 



to as Vaughan-Preston (VP) gap. lHartmann et al 
( 1984 ) explained it as a combined effect of CA satura- 
tion in active stars and a basal level in the CA indicator 
used by VP, due to the photospheric flux. These, it was 
claimed, enhance the impression of a gap. A variation 
in the local stellar birthrate was also invoked. Other 
studies consider the VP gap a result of the nature of 
the C A evolution (e.g. IDurnev et al.lll981t iMiddelkoop 
1982). Our Figure [1] corroborates the latter hypothesis: 



all stars younger than 1.2 Gyr, except one, lie above 
the gap, all stars older than 1.4 Gyr with the exception 
of one, lie below it, indicating a drop of CA level in a 
very short time. The two exceptions are probably due 
to a particular phase of the activity cycle. It is worth 
noticing that the stars in either side of the gap differ 
not only in CA level, but also in its trend as a func- 
tion of temperature: logi?^ (or equivalently \ogR' HK ) 
depends weakly on the temperature for stars above the 
gap, while it has a distinct decreasing trend for stars 
below. This can be seen from our Figure [TJ and it ap- 
pears clearer when a larger t emperature range is con- 
sidere d, like, for instance, in IrVlamaiek fe Hillenbrand 



( 200cj see Figure 4 therein). Furthermore, short-term 
temporal variations of CA are large and irregular for 
active stars an d small and regular for inactive ones 
( Vaughan 1980h . There must be a major event at a 
given time of the stellar main sequence life time, that 
changes the way radiative heating mechanisms occur in 
the chromosphere, and, as a consequence, their depen- 
dence on stellar parameters as well as the shape and 
le ngth of CA cycles. 



Fawzv et al.l (|2002f) present theoretical calculations 



that reproduce the observed trend of CA with stel- 
lar temperature. They invoke two heating mecha- 
nisms: the magnetic-wave and the acoustic-wave heat- 
ing. Theoretical chromospheric fluxes for the former 
mechanism match the observed flux of the active stars, 
while theoretical fluxes for the latter mechanism match 



t that stars occupy two very distinct sequences in a 
rotation period-versus-cycle period diagram. Not all of 
these works can be used to explain our results. What is 
relevant for the present discussion is the possibility of a 
change in the nature of the dynamo mechanism taking 
place at a specific point of the main-sequence lifetime 
of solar-type stars. 

The most important conclusion of our analysis is 
that, before and after the fast drop of CA from above 
to below the VP gap in only 200 Myr, its alleged 
smooth decay must be much less imp ortant than short 
term variations. We noted first in iPace fc Pasauini 
( 20041 ) that CA in intermediate age clusters had al- 
ready dropped to the solar level. After the reanaly- 
sis of the data using temperature determined spectro- 
scopically from iron lines, 1 out of the 7 stars in the 
two intermediate age clusters IC 4651 and NGC 3680 
turned out to lie in the VP gap. In addition, the 
membership of one of the two stars in NGC 3 680 has 
been questioned Anthony- Twarog et al. ( 20091 ). How- 
ever, the data still point undoubtedly to a lack of evo- 
lution after 1.4 Gyr, since 5 out of 6 secure interme- 
diate age stars have CA levels equal or lower to those 
spanned by the Sun and by M 67 stars (empty dots 
in Figure [TJ This conclusion is al s o corr oborated by 
the work of iLvra fc Porto de Meliol |2005l) . The other 
part of the conclusions, that regarding the time in- 
terval between 0.7 and 1.2 Gyr (filled dots in Fig- 
ure [JJ, was achieved after the analysis of 5 stars in 
NGC 5822 and IC 4756, together with the older data 
on Hyades and Praesepe ([Pace et al.ll2009l ). We could 
show that metallicity is unlikely to play a major role, 
since there is no significant difference between three 
almost coeval open cluster s at different me tallicitics: 
Praesepe, at [Fe /H]=0.27 dPace et alJbOQSl ). Hyades, 
at [Fe /HI =0. 13 (IPaulsori et al.ll2003l ). and IC 4756 at 
[Fe/H]=0.01 |Pace et al.ll2010l) . The quoted iron abun- 
dance for Praesepe is s omewhat different from othe r 
published works (see e.g. Anthony- Twarog et al. 20091) . 
However, this result is based on the high-quality spec- 
tra described in Section [JJ the iron abundance was 
obtained adopting both spectroscopic and photometric 
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Hyades, 0.7 Cyr 
Praesepe, 0.7Gyr 
IC4756, 0.8Cyr 
NCC5822, 1.2Cyr 
Stars in the gap 
NGC3680, 1.4Cyr 
IC4651, 1.7Cyr 
M67, 4.3Cyr 
Sun, 4.5 Cyr 



Fig. 1 Diagram of CA as a function of the temperature using log R' K as a proxy for CA . For the Sun, 2 datapoints, 
corresponding to the maximum and the minimum of an activity cycle, are plotted. Datapoints relative to stars above 
Vaughan-Preston gap are plotted with empty symbols, datapoints below with full symbol s. The typical la er ror is 0.1 dex 
on the ordinate and 110 K on the abscissa. The ages indicated in the legend are taken bv lSalaris et al.l (| 20041 ). 



temperatures, several ways of normalizing for the con- 
tinuum were tried, and even different people performed 
EW measurements in order to avoid personal bias, and 
in any case the result was [Fe/H] significantly higher 
than 0.2 dex. 

Our conclusions on CA evolution are still based on 
few sta rs per cluster and they partly rely on the age 
scale of Salaris et al. ( 20041 ). which uses a calibration of 
a photosphcric age indicator made on 11 clusters. In 
order to strengthen our results, or to disprove them, 
we need more CA measurements in solar-type stars in 
open clusters based on high-resolution spectra and di- 
rect age determinations of the open clusters used based 
on improved colour-magnitude diagrams. We should 
both reobserve the already analysed stars, to obtain a 
time-averaged CA level, and choose other targets. 

The scenario suggested here contradicts the common 
belief that CA is well correlated with age, however our 
investigation does not include stars younger than the 
Hyades, some of which are indeed more active than the 
Hyades stars. In addition, the fact that young stars 
are active and old stars are not is not questioned here, 
and it causes a weak correlation between CA and age. 
We believe that data used in the literature to prove 
a deterministic relation between CA and age do not 



disprove our alternative view, and that more data are 
necessary to ac hieve a definitive conclus ion. 

For instance ISoderblom et al.l (|199lf ) used Mount- 
Wilson data about several solar-type stars in visual bi- 
naries and single F dwarfs, with ages from Stromgren 
photometry in order to prove that the relationship be- 
tween CA and age is, using their words, determinis- 
tic and not statistical. Figure 3 in the aforementioned 
paper shows that log R' HK and the logarithm of age 
are indeed correlated. However we notice, analysing 
data in Table 1 and 2 therein, that among stars with 
\ogR' HK > —4.9, i.e. more active than the Sun, the 
correlation between \ogR' HK and age is not significant, 
the Pearson correlation coefficient is about 0.1. As for 
the other stars, i.e. those with logi?^ A - < —4.9, the 
Pearson coefficient is 0.35, and its sign is positive, un- 
like what is expected according to an activity decay. 
In either group of stars there is, instead, a significant 
correlation between log R' H K and the B-V colour, and 
it has opposite signs. 

Ano ther calibration o f CA evolution with time is 
that of Lachaume et al. (|l999h . From Fig. 4 therein, 
it can be seen that their result agrees with ours as 
far as inactive stars are concerned, i.e. CA does not 
evolve after it has crossed the VP gap. As far as active 
stars older than the Hyades are concerned, there are 
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4 to 6 data-points (for two stars it is not possible to 
say whether or not they are younger than the Hyades 
and therefore out of the range considered in the present 
study). For this group, the Pearson coefficient indicates 
a level of anticorrelation between age and log R' H K that 
is weak (-0.27) or fair (-0.57), depending on how many 
stars we consider. 



From iMamaiek fc Hillenbrand! (|2008[ ) it is clear that 
the one open cluster that supports a strict monotonicity 
of CA time evolution (once short-time scale variations 
are smoothed out) in the age range between that of the 
Hyades and that of M 67 and the Sun, is N GC 752. This 
cluster, according to ISalaris et al. ( 2004 ). is older than 
NGC 5822 and younger than NGC 3680, i.e. exactly 
in the range i n which we expect the tr ansition to oc- 
cur. However, IMamaiek k, Hillenbrand! also show that 
companions of a binary system tend to have similar CA 
level, especially if they have similar colours. According 
to the assumption that CA is a good age indicator, this 
circumstance is easily interpreted as due to the fact that 
companions of a binary system are coeval. We suggest 
that it could be instead related to the nature of the 
binary systems. 

We conclude that data in Soderblom et al. ( 199lh: 
Lachaume et all |l999[ ); IMamaiek fc Hillenbrand! (l2008l )| 



are compatible with the conclusion that, within the age 
range from 0.7 to 1.2 Gyr and from 1.4 Gyr to solar age, 
any age-activity relationship must be weaker than the 
short term variations 
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